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Abstract 

We study the discrete massless Gaussian Free Field on Z"*, d > 2, in presence of two types 
^ of random environments : (1) (5-pinning at height of inhomogenous i.i.d. BernouUi strengths; 

(2) square-well potential supported on a finite strip with i.i.d. Bernoulli reward/penalty coeffi- 
I cients e. We prove that the quenched free energy associated to these models exists in M+, is 

^ self-averaging, and strictly smaller than the annealed free energy (whenever the latter is strictly 

positive). Moreover, for model(2), we prove that in the plane (Var(e), E(e)), the quenched critical 
line (separating the phases of positive and zero free energy) lies strictly below the line E(e) ~ 0, 
showing in particular that there exists a non trivial region where the field is localized though 
I— I repulsed on average by the environment. 

(-H Keywords : Random interfaces, random surfaces, pinning, disordered systems, Gaussian free field. 
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1 Introduction 
1.1 The models 

We study two models belonging to the so-called effective interface class. For A being a finite subset of 
Z"^ (denoted by A (s Z"^), let if = {fx)xeA represent the heights of the interface above or below sites in 
A. The interface is thus living in Z'^^^. 




Figure 1: A 2-dimensional interface. 

The variables ipx can also be seen as continuous unbounded (spin) variables, we will refer to ip as "the 
interface" or "the field". The models are defined in terms of a Gaussian pair potential and a random 
potential term at interface height zero (or close to zero). 

In the first model, we consider only attractive potential. The so-called (5-pinning measure is given by : 
Ulidif) = ( "^{"Px - ^y? j n (^'^^ + '/^(^xSoid'fx)) n '^o(d¥'y), (1) 

^ \ xr^y ) a:eA J/eA'= 

where /3 > is the inverse temperature, A d Z"', x ~ y denotes nearest neighbour relation in A U 9A 
and (5o is the Dirac mass at height 0. e := {ex^xf^K is given by positive i.i.d. Bernoulli random variables 
on Z"^ standing for the quenched disorder, which thus represents a random attraction strength at height 
0. Z\ is the partition function, i.e. it normalizes [i\ so it is a probability measure. 
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In the second model, we consider attractive/repulsive potential, and define the following measure 



Pl,a{df) = exp I ^i}Px - ^yf + exl[^^(.[^aM] I n ^"^^ n ^o{d^y), (2) 



where Ij] denotes the indicator function, a > and e is again a Bernoulli product random field on Z , 
which is now allowed to take negative values, thus repulsing the interface. 

Note that both models contain two levels of randomness. The first one is e which we will refer to as 
"the environment". The second one is the actual interface model, depending on e. By performing a 
standard scaling we can eliminate the inverse temperature /3, thus henceforth we assume /3 = 1; cf. 
Remark O 

The dimensions 1 and 2 are physically relevant as interface models. In this paper we focus on d > 2 



since 1-dimensional models have been well-studied in the last decade (see Section (1.3.1) below). 
The questions we are addressing in this framework are the usual ones concerning statistical mechanics 
models in random environment : Is the quenched free energy non-random ? Does it differ from the 
annealed one ? Can we give a physical meaning to the strict positivity (resp. vanishing) of the free 
energy ? What can be said concerning the quenched and annealed critical lines (surfaces) in the space 
of the relevant parameters of the system ? 

For models defined by ([T]) and (|2]), we prove that the quenched free energy exists and is non-random. 
It is always strictly positive in case of model ([T| , and strictly smaller than the annealed free energy for 
both models. 

For model (|2]), we also investigate the corresponding phase diagram : in the plane (Vare(e), E(e)) (i.e. 
the variance and expectation of e), we prove that the quenched critical line (separating the phases of 
positive and zero free energy) lies strictly below the line IE(e) = 0. Thus there exists a non trivial 
region where the field is localized though repulsed on average by the environment. 
We do not treat here the questions concerning critical exponents of the system, nor the order of 
the phase transition, in presence (absence) of disorder. This has been done for a certain class of 
1-dimensional systems (see below), and is a much more difficult problem in dimension 2 and above. 

1.2 Known results about homogenous models 

Our work relies on |5],j9]. They focus on the case Cx — const which we will refer to as the homogenous 
pinning model. A review article about localization and delocalization of random interfaces (in a non- 
random environment) can be found in [23j. We briefly describe some important results here. 
Let us now consider the model ([l]) [^with Cx = £ with e > for all x G Z*^, and denote the associated 
measure by /x^*^, the superscript standing for boundary condition. The so-called free measure is 



0,0 /^0,0\-l 

/^A = (^A ) exp 



\ xr^y J xeA j/eA= 



is centered Gaussian, and thus characterized by its covariance matrix. There exists a useful represen- 
tation of the latter in terms of the Green function of the simple symmetric random walk on Z'^ 



^l]f{ipx^y) = Ex 



■ta-I 
. n=0 



(3) 



^Similar results hold for the homogenous version of model (^2| with = e > Q for all x £ Z*^, see [5] 
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where, {{Xn)n, Px) is a simple symmetric random walk on Z'^, starting at x, and t\ is its exit time 
from A. As the simple random walk is recurrent in d = 2 and transient in d > 3, the variances diverge 
in d = 2, and stay bounded for d > 3 as A f Z'*. We refer to |21| for random walks estimates. 
To analyze the properties of the model with pinning (e > 0), it is convenient to map it onto a model of 
a random walk in an annealed random environment of killing obstacles. To achieve that, we can expand 
the product H^eA + ^^oidfx), as in a high temperature expansion. Let A = {x € A : (px = 0} 



be the random variable describing the set of pinned sites. By simple calculations (see (16) for more 
general case), for any measurable function / : i— )• M, we have : 

/^f (/) = E ^'f(f I = ^) • = ^) = E ^a\a(/) • = A). (4) 

AcA AcA 

For simplicity we denote i^%{A) := ^^^^{A = A). Conditionally on A, /U^*^ is a Gaussian measure with 
covariances described by the random walk representation above. 

Important results concerning the pinned sites distribution were obtained in |5|. First of all, the measure 
is strong FKG in the sense of [TO]. Moreover, it can be compared with i.i.d. Bernoulli distributions. 
Let be the Bernoulli product measure (or site percolation) of parameter a £ [0, 1] in the box A, 
namely the measure on subsets of A given by : 

B^(^ = ^) = al^l(l-a)l^l-l^l 

By [5l Theorem 2.4] there exist constants < c_(d) < C-^-{d) < oo such that for any A, any B C A, 
and e sufficiently small in d = 2, we have 

^ 0) < u%iAnB = 0) <M''+^'^^^^'\AnB = 0), (5) 

where 

, , r e|loge|-i/2 d = 2, 
'^'^ = [e d>3. 

For d > 3, an even stronger statement is true : is strongly stochastically dominated by B^^'^^'' and 
strongly stochastically dominates B^"*-*^^^, i.e. for any C C A: 

c-{d)e < i^iix G A I A\{x} = C)< c+{d)e. (6) 

Concerning the behavior of the interface, it is known that an arbitrarily weak pinning e is sufficient 
to localize the interface. Indeed, in [9], Deuschel and Velenik proved, for a class of models including 
ours, that the infinite volume Gibbs measure /z^'" exists in all d > 1 and that for any e small enough 
and all K large enough 

( K d = 1, 

-\ogif^\^o> K)^A KyiogK d = 2, (7) 

( d>3. 

The so-called mass, or rate of exponential decay of the two-point function, associated to the infinite 
volume Gibbs measure /x^'" is defined, for any x G S'^^'^, by 

m%x) := - lim jJogCoVe{<fo^[kx])- 



a Xd b means that there exist two constants < ci < C2 < oo, depending only on d, such that cib < a < C2b. 
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where [x] is the vector of integer parts of x's coordinates. In fTTf loffe and Velenik showed that for any 
d>l, 

inf m^x) > 0. (8) 

The locahzation of the interface becomes weaker as e — )• 0, we can quantify this by studying the 
behavior of the variance and the mass of the field in this hmit. The most precise results were proved 
by Bolthausen and Velenik in |5], and can be stated as follows: 
For d = 2 and e small enough, 

= ^\loge\ + 0{log\loge\) (9) 
vr 



For e small enough, 



1.3 Known results about disordered models 
1.3.1 Models on Z 

In [2], Alexander and Sidoravicius studied the 1-dimensional (attractive/ repulsive) model. They 
consider a polymer, with monomer locations modeled by the trajectory of a Markov chain in 
the presence of a potential (usually called a "defect line") that interacts with the polymer when it visits 
0. Formally, let Vi be an i.i.d. sequence of 0-mean random variables, the model is given by weighting 
the realisation of the chain with the Boltzmann term 

exp (^^Yl{u + Vi)lix,=o^ ■ 

Their purpose was to study the localization transition in this model. If a positive fraction of monomers 
is at 0, we say that the polymer is pinned. In the plane /3 vs. u, critical lines are defined : for 
(3 fixed, let Mc(/3) (resp. u^(/3)) the quenched (resp. annealed) critical value of u above which the 
polymer is pinned with probability 1 (for the quenched (resp. annealed) measure). They show that the 
quenched free energy and critical point are nonrandom, calculated the critical point for a deterministic 
interaction (i.e. Vi = 0) and proved that the critical point in the quenched case is strictly smaller. 
Note that when the underlying chain is a symmetric simple random walk on Z, the deterministic critical 
point is known to be 0, so having the quenched critical point (uc{P)) strictly negative means that, even 
when the disorder is repulsive on average, the chain is pinned. This result was obtained by Galluccio 
and Graber in for a periodic potential, which is frequently used in the physics literature as a "toy 
model" for random environment. 

In |16| . Giacomin and Toninelli investigated the order of the localization transition in general models 
of directed polymers pinned on a defect line. 

They prove that for quite a general class of models, as soon as disorder is present, the transition is 
at least of the second order, i.e. the free energy is differentiable at the critical line and the order 
parameter (which is the density of pinned sites) vanishes continuously at the transition. 
This is particularly interesting as there are examples of non-disordered systems with first order transi- 
tion (cf. for example Proposition 1.6, for (1 -|- d)-dimensional directed polymers and d > 5). The 
result thus implies that the introduction of a disorder may have a smoothening effect on the transition. 
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For 1-dimensional models, the renewal structure of the return times to plays important role, in par- 
ticular it simplifies a lot of calculations. In Alexander emphasized this fact by assuming that the 
tails of the excursion length between consecutive returns of X to are as n~^(j){n) (for some 1 < c < 2 
and slowly varying 0). He analysed the quenched and annealed critical curves in the plane {u, (3) for 
different values of c, showing that for c > 3/2 at high temperature the quenched and annealed curves 
differ significantly only in a very small neighborhood of the critical point, whereas for c < 3/2 the 
quenched and annealed critical points are equal. This was a prediction made by theoretical physicists 
on the basis of the so-called Harris criterion (see [13], Section 5.5, for more informations). The relevant 
case in the framework of this paper is the case of the Markov chain given by a simple symetric random 
walk on Z, which corresponds to c = 3/2 and </>(n) ~ K for K > 0, which is borderline. The author 
is unable to say whether the critical behavior is altered by the disorder in this case, but progress on 
this question has been made recently by Giacomin, Lacoin and Toninelli in |14| and |15) . They prove 
that in the borderline case c = 3/2, the disorder is relevant in the sense that the quenched critical 
point in shifted with respect to the annealed one. They consider i.i.d. Gaussian disorder in the first 
paper and extend the result to more general i.i.d. laws, as well as refine the lower bound on the shift, 
in the second paper. Note that this case includes pinning of a directed polymer in dimension (1+1) 
as already mentionnecQ but also the classical models of two-dimensional wetting of a rough substrate, 
pinning of directed polymers on a defect line in dimension (3 + 1) and pinning of an heteropolymer 
by a point potential in three-dimensional space. 

1.3.2 Models on U^, d>2 

The only result about random pinning models we are aware of is [18] . In this paper, Janvresse, De La 
Rue and Velenik considered the model (jl) in dimension 1 and 2 where e € {O,??}^'', which models a 
interface interacting with an attractive diluted potential. They show that the interface is localized in 
a sufficiently large but finite box (in the sense that there is a density of pinned sites) if and only if 
the sites at which the pinning potential is non-zero have positive density. Note that in this paper they 
characterize the set of realizations of the environment for which pinning holds (the disorder is fixed, 
not sampled from some given distribution), which is stronger than an almost sure result. 
We also mention a series of papers by Kiilske et al. (|19). |22| . |20) which study a model with disordered 
magnetic field (instead of disordered pinning potential ). For example, Kiilske and Orlandi studied the 
following model in dimension 2 



e.O,(»?)/T N 1 



V x^y xgA / xgA y&A" 

where {rix)x£A is an arbitrary fixed configuration of external fields and V is not growing too slowly at 
infinity. Without disorder (?] = 0), the interface is localized for any e > One could expect that 
in presence of disorder and at least for very large e the interface is pinned. However, the authors show 
that this is not the case : the interface diverges regardless of the pinning strength. This implies that 
an infinite-volume Gibbs measure for this model does not exist. One could hope for the existence of 
the so-called gradient Gibbs measure (Gibbs distributions of the increments of the interface). In |22] 
Van Enter and Kiilske proved that such (infinite volume) measures do not exist in the random field 
model in dimension 2. Note that gradient Gibbs measures may exist, even when the corresponding 
Gibbs measure does not. This happens when the interface is locally smooth, although at large scales 
its fluctuations diverge. This is the case for the two-dimensional Gaussian free field. 



We speak about "polymer in dimension (d + 1)" when the state space of the Markov chain X is of dimension d. 
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1.4 Open problems 



The model studied in this paper lends itself to number of extensions. We list here a selection, with 
brief comments. 

• Path-wise description of the interface. In the positive free energy region, for both of 
our models, one expects localization, i.e. the finite variance of and exponential decay of 
correlations. A much more difficult question concerns the behavior of the interface near the 
critical line. Does it behave the same as in the homogenous case (i.e. second order transition 
with the density of pinned sites decreasing linearly for (i > 3, and with a logarithmic correction 
for d = 2) ? Or does the presence of disorder have a smoothening effect on the transition (as it 
was proven for certain 1-dimensional models) ? 

In the zero free energy region (in the attractive/repulsive setting), we expect the behavior similar 
to the entropic repulsion for the GFF : in a box of size n the interface should be repelled at height 
lb log n in d = 2 and zt-y/log n in d > 3 (see |8j and references therein for details). The it stems 
from the fact that our model is symmetric with respect to reflection at zero height, hence with 
probability 1/2 it either goes upwards or downwards. 

• Description for non Gaussian case. A natural conjecture is that the behavior of the model 
is the same if we change the Gaussian term (ipx — ^y)"^ to any other uniformly convex potential 
V{ipx — ^y)- Let us note that Griffiths inequality is not proven even in the homogenous case, 
which can make the problem difficult to handle. 

• Non-local interactions. We restricted our work to the case of the nearest neighbors inter- 
actions. We suspect that the results holds true for fast decaying interactions, at least with 
condition like in [5, (2.1)] (which ensures a control of the random walk's behavior in the random 
walk representations). As the behavior of the homogenous pinning model beyond this regime is 
not known, we are unable to pose any further conjectures. 

• Non i.i.d. pinning laws. Going beyond the i.i.d. case is a very interesting direction. Two 
natural cases would be the stationary Bernoulli field or the chessboard like configuration. These 
questions may be closely connected to convexity /concavity properties of the free energy function 
in the homogenous case. The understanding of this case is still limited. It would be interesting 
to know if finite range environment laws change the picture, as it is sometimes the case in models 
with bulk disorder but it seems difficult to answer to this question rigorously. 

• Geometry of pinned sites. The geometry of the pinned sites is still not fully understood in 
the homogenous case. For the d> 3 the law of pinned sites resemble a Bernoulli point process. It 
is conjectured that once the pinning tends to zero, under suitable re-scaling, this field converges 
to Poisson point process. For d = 2 the situation is not clear at all, though it is expected that 
the dependency between the points will be preserved in the limit (implying the limit being non- 
Poissonian) . 

Not only these questions propagate to the non-homogenous case but also new ones arise. E.g. for 
the attracting/repulsive model it would be interesting to study the joint geometry of attractive 
and repulsive sites. 

• Models with wetting transition. The effects of introducing a disorder in other models with 
pinning might be interesting, for example in models exibiting a wetting phenomenon. In the case 
of the massless Gaussian model in d = 2, it is known [6] that the wetting transition takes place 
at a non-trivial point. A natural question to ask is, if adding disorder shifts the transitions point. 



7 



2 Attractive environment 



2.1 Introduction 

2.1.1 Description of the results 

We denote by Q := M^'* the set of configurations of tlie model. Let A be a finite subset of Z'^. Tlie 
finite volume Gibbs measure in A for the discrete Gaussian Free Field with random 5-pinning e, with 
boundary condition, is the probability measure on defined by : 

/^fid^) = ^e-^'"'^^'"^ n (^'^^ + V^^-^oidiPx)) n ^oid^y)^ (12) 
A xeA yeA<= 

where A'^ := Z'^\A and 6o is the point mass at 0. The normalization constant is the so-called 
partition function of the model. The Hamiltonian in A is given by 

^aM = ^ E (^--V>yf, (13) 

{x,y}nAj^0 
x^y 

where x ~ y means that x and y are neighbouring vertices in U^. The non-homogenous pinning 
environment is denoted by e := {ex)x<^A- We consider the law of e given by a Bernoulli product 
measure. Namely, independently on each site x G A : 

e with probability p , . 

e with probability (1 — p), 

for some fixed triplet {e,e,p) such that < e < e and p G (0,1). In this section we consider only 
positive values of the environment, meaning purely attractive environment. Let 

e* := pe + {l- p)e, and (7^ := pe^ + (1 - p)e^ - (e*)^. (15) 

be the average and variance of e^, respectively. 



Remark 2.1. 1. The measure (12) is the weak limit of the measures /x^^n, given by once we 
choose (a")n>i, (e")„>i such that a" — )• and 2a"(e^" — 1) = e^. 

2. The temperature parameter enters only in a trivial way. If we replace the field {(px)xGA by 
{■\/P(f)x)x£A, o,nd {ex)x<^A by {\^ex)xeA we have transformed the model to temperature parameter 
/3 = 1. Thus in the sequel, we will assume (3 = 1. 

By definition the annealed model is obtained by averaging the Gibbs weights over the disorder; the 
annealed model corresponding to environment e is thus the same as the homogenous model with 
pinning strength e*. 

We define the quenched (resp. annealed) free energy per site in An := {0, ...,n — l}*^ by (e) (resp. 

flie) := n'^log 1^^^ , /^(e) := n-'^log {^^^ . 

where denotes the partition function of the model with no pinning (i.e. e = e = 0). This 

normalization is chosen such that the free energy of the free model is zero. 
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In the next section we show that (/^ converges to a non-random quantity: the infinite volume 

quenched free energy, denoted by /''(e), i.e. 

/^(e) ^ /'1(e) = Ee(/'i(e)) a.s. and in as A t provided |aA|/|A| ^ 0. 

The same proof can be used to show that /^(e) — )• /'^(e) as A t Z"^ provided |5A|/|A| — )• (see [T6] 
for the proof in d = 1). Comparing the free energy for the quenched and annealed models is useful 
for studying the effect of the environment on the system. The case when they are different indicates 
that the disorder of the environment has some macroscopic effect. By the Jensen inequality (applied 
to log) and convergence of /^(e), it is easy to present that 

ne) < ne). 

We recall that in this section we consider purely attracting environment. Thus we have /''(e) > for 



all {e,e,p) such that < e < e and p G (0, 1), as we prove in Fact 2.5 Moreover, the density of pinned 
sites in the box A„, is strictly positive as n tends to infinity, as we prove in Fact |2.6| In Theorem |2.8[ 
we show the main result of this section, namely that /''(e) < /'^(e) for all {e,e,p) such that < e < e 
and p G (0, 1). We provide also some estimates on the difference f^{e) — /''(e) in terms of = Var(e). 



2.1.2 Extension of pinned sites representation to inhomogenous pinning 

In this section we extend the decomposition Q to our inhomogenous model. We write 
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A AcA \x^A ?/eA\A / JieA': 

Efn«.)}IO'i^). (16) 

AcA VxeA ' 



= :vl{A) 

In other words v\ is the marginal of the measure \x'^(dip) giving the distribution of pinned points. We 
can obtain a formula for the ratio of partition functions appearing in the free energy : 

^A AcA VxeA / ^A 

2.2 Existence of quenched free energy 

Theorem 2.2. For d > 2, let /^(e) := |A|-Mog(Z2' /^°'°) he the quenched free energy per site in 
A (S Z'^. Then, the limit 

/"(e) := lim /^(e), 

AtZ'* 

exists almost surely and in L?' , and does not depend on the sequence A f Z'^ provided it satisfies 
|5A|/|A| — )• 0. Moreover, /''(e) is non-random, i.e. 

/''(e)=Ee(/''(e)) a.s. 
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Proof. We prove the existence of the hmit along the sequence of boxes Bn '■= A2"-i. The generalization 
to all sequences A f Z'^ such that |9A|/|A| — )• is rather standard (cf. for example |24)). 
Let us write aefl{dip) := nxe_B„(^Va; + ^x^oidifx)) Y[y^B„ ^o{d<-Py)- We recall also notation ( [l3| ) for the 
Hamiltonian. We will cut Bn in 2°^ sub-boxes denoted by sjl^. Let X := (U-=i dB^^l^)\dBn be the 
interface between the sub-boxes. In order to prove the existence of the limit along {Bn)n, we first derive 
a "decoupling property". Namely there exist c„ > such that c„ < oo and l'^^'^^ — 11^=1 'il — '-'^ 

for any realisation of e, where e^*) is the restriction of e to the box B^}_^. This allows us to prove that 
expectation of /^^(e) converges, and its variance tends to zero, i.e. /^^^(e) — )• c € M almost surely, in 
and in L^. 



Lower bound on 

We have. 



where e^*-* denotes e restricted to -6^*1 1. Hence, 



for some constant Cmin > 



n (e) = 2-"'^ 



Upper bound on 

We will first prove that with high probability < 2"*^ for alH G X and some small (5 > 0. This will 
allow us to "force" if to be zero on X for small energetic cost. 

Lemma 2.3. There exists Ci, C2 > such that for all i £ An and for all n sufficiently large, 

>T) <Cle-^^^'/l°g^ 



Proof. Using (16), 

^^:'^{\^.\ >T) = ^n(^) • iV^i^^i > ^) 

ACAn 

Now, fJ'^jf\^j^ is Gaussian, therefore, ipi ~ AA(0, cta) with 

aA := Var°'°^^(99i) < Var°'°(99o) < Clogn 
for some C > and n large (cf. [7]). We can use the Gaussian tail estimate : 

In the sequel, the notation C, C , C" will be used for positive constants that may change from line to 
line. L et us define X := X [J {dX n -Bn), which is a thickening of X consisting of 3 "layers". Then, 
Lemma 



2.3 



allow us to control the height of the field on X C Bn '■ 
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Hence, 



Jifi\j^el~2^",2^"] 



v^-C2'5"^/-^//„2"(2■5+<^-l)^5n2"(<^-l) /7e(*),0 

1=1 

This leads to 



2*^ 



2d 

1=1 



Combining our two bounds, we obtain : 

C^in2-" < f^Je) - ^l;/l_,(e«) < C^ax2"(2^-i) (19) 



1 = 1 

for some constants Cmim C'max > 0. 

Expectation of (e) converges and its variance tends to zero 

By (19), it is easy to see that : 

1 

|Ee(/^Je)) -Ee(/^_^(e))| < |Ee(/^Je)) - Ee(^ /^_^ (e«))| < 02^^''-'^ (20) 

i=l 

The right hand side is summable for 6 < 1/2 hence Ee(/^ (e)) converges as n — t- oo. Using indepen- 

(i) 

dence of environment among the boxes we can deal with the variance. Let us write 

i=l 

where if„ is the above error term, \£n\ < C'2"(2<5-i). Then, 

2d 2'' 

Vare(/^^ (e)) = Vare(2-'^ ^ /^_^ (e^)) + Vare(^:„) + 2CoVe(2-'^ ^ /^_^ (e«), £:„) 

1=1 i=l 

= 2-'^Vare(/^_^(e)) + Vare(^:„) + 2CoVe(/^_^ (e), £:„) 
< 2-'^Vare(/^_^(e)) + Vare(f„) + C"2"(2^-i)Ee(/^_^ (e)) 
Now, since (Ee(/^ ^(e)))„ converges, we get the following upper bound on the variance: 

Vare(/^Je)) < 2''^Vare(/l^_^ (e))+C22"(25-i) + C''2"(25-i) (21) 

We deduce Vare(/^^(e)) 0, and conclude /^^(e) c G M in L^. The sequence converges also almost 
surely, since all error terms are summable. □ 
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Exactly the same proof until ( 19 ) yields 



Theorem 2.4. For d>2, let /^(e) := |A|-i log(EeZ^' /Z"'") be the annealed free energy per site in 



A d Z'^. Then, the limit 



r(e) := lim /^(e) 



exists and does not depend on the sequence provided it satisfies |5A|/|A| — )■ 0. 

2.3 Bounds on the quenched free energy 

2.3.1 Strict positivity of the quenched free energy 

In this section we prove that the free energy of the model is strictly positive for all environments 
(e, e,p) such that < e < e and p G (0, 1). That implies the localization of the field, namely the strict 
positivity of the limiting density of pinned sites. Our results are valid for all > 2. 

Fact 2.5. For all (e, e, p) such that < e < e and p G (0,1), we have /^(e) > 0. 



Proof. Using (17), and the fact that < e < e, we can write 



r(e) = „lim n-'^log ( J] J] -^^] ^ il.- -"'1°^ f E = 

\AcAnxeA ^An / \AcA„ ^An / 

Moreover, it follows from [5] Theorem 2.4, that the asymptotic density of pinned points under fj,j^ is 
stictly positive uniformly in A„ for n large enough. This immediately implies by definition of the free 
energy that /'^(e) > 0. □ 

Fact 2.6. Let < e < e and p G (0, 1) then there exists c > such that for almost every realization e 
we have 



4v.=o] >cj =1. 

Remark 2.7. More is known for d = 1. By [2] there exists p = p{e*, a) such that 

L[^,=0] e{p-e,p + e)\ =1. 



lim lim u.' ( n Ir 

\ xeAn 

Proof. Let S > 0, and As := {YlxeAn '^['Px=o] ^ Sn'^}. Let us fix some e > 0. By definition we have 



^A„ ^A„ xeA 



Using the pinned sites decomposition (16) we get : 



"^A„ AcA„ ^A„ AcA„ ^A„ 



where e^ := W^^^a^x, 1 = (7x)xeA„ and 7^; := e^/e. Similarly, 



xeA 
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We thus have 



n-'' log fifjAs) < S\og{e/e) + n'^^log - n-'^log 

^^n Jin 



By Fact 2.5 the last term converges to a strictly positive quantity. Moreover, it is not difficult to show 
that < + Ce)" for some C > uniformly in n, hence the second term tends to uniformly 

in n as e — )■ 0. Choosing e > and 5 > small enough we get 

limsupn^'^log^^''^(A5) < 0. 
and the result follows. □ 

2.3.2 Strict inequality between quenched and annealed free energies 

Theorem 2.8. Take (p,e, e) such that p £ (0,1) and < e < e. Let e be an i.i.d. Bernoulli envi- 



ronment of parameters {p,e,e) as described in (14)- We recall the notations e* = pe + (1 — p)e and 
0-2 = e2p + 6^(1 -p) - e*^. 

(a) Letd>3. Then, 

Pie) < r(e). 
Moreover, there exists c^{d^p) > such that, 

f-{e)- r{e)>C3id,p)a + 0{a^) as a ^ 0. 

(b) Let d= 2. There exists e > uniform in p,e,e such that if e* < e, then 

r(e) < r(e). 

Moreover, under the same hypothesis on e* , there exist C2{p,e*) > such that, 
r(e) - /'1(e) > C2{p, e*)a + 0{a^) as a ^ 0. 



Proof. By Theorem 2.2 we have 

/ 76,0^ 



/q(e)-/-(e) = lim n-'^Eelog 



n— >oo 



By ( 22 1 we conclude that our goal is to prove 

lim sup n'^'lEe log [y''^ (7'^)) < 0, 

where 7^. := Cx/e*. 

We will now outline the proof strategy. We notice that if was a Bernoulli product measure with 
some intensity A G (0, 1) then the theorem would follow easily. Indeed, in such a case, we can calculate 
directly the expectation under log and the expression in question is equal to 

n-'^Ee log I n (^^^^ + 1 - ) = log(p7o + 1 - p) < 0. 

\a;eA„ / 
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The last inequality follows by the strict concavity of the logarithm and the Jensen inequality. Our aim 
is to use ^ in order to make this argument rigorous. 

Let us now introduce an additional randomization. We will use the uniformly distributed permutations 
of An, writing the corresponding expectation by E. It is easy to check that for any i.i.d. pinning law 
EeE(-) = Ee(-). Let TT be a permutation of the n'^ vertices of A„, chosen uniformly at random. By the 
Jensen inequality, we have 



n-'^Eelog(z.^*(7-^)) = n-'^EeElog^.;^M Ht.W 



< n-'^Ee log Ez^f ( J] 7,(,) ] = n-^Ee log Ez^f (7-(-^) 
\ieA J 

where it{A) := {vr(z) : i £ A}. Intuitively, Ei/^ is the expectation of the distribution of pinned sites 
"scattered" by a random permutation. Thanks to this we can work with a uniformly distributed set 
of pinned points, provided we know its cardinality. More precisely, 

k=0 ^ ^ \AcA„:m=fc / 

Our aim is now to use stochastic majoration/minoration by Bernoulli measures to obtain an upper- 
bound for (1^1 = A;). Difficulties will follow from the fact that {|^| = A;} is not an increasing 
event. 

Case (i > 3. We recall the notation and the stochastic dominance results introduced in ([6|. In d > 3, 
we know that stochastically dominates a Bernoulli product measure with A := C-{d)e* > 0. 
For any set A C A, we have 

B^(^ = ^) = aI^|(i-a)"'"I^i 

Let 6n,A:(A) := ~ A)"''"'^. As the event {|^| < k} is decreasing, we have the following upper- 

bound : 

<'{\A\ = k) < .^;*(|.4|<t)<Bj(|.4|<*,-) (23) 



. . ^K,x{k) ' \ ^ -"--L 6„;^(i -M) 



(*) 

Now, for i < [An^J , 



bn,x{i) _ i + l < 1 



bn,x{i + 1) Wi-i X 
This gives an upper bound for k < [An'^J , namely 

i^U\A=k)<cin''-M^ni\A\=k). 

For k G [ [An'^] , n'^] , define Afc := kn^'^. As by the Stirling formula bn^x^ik) > C2n~'^/'^ > uniformly in 
k, we can use the bound: 

u^'{\A\ = k)<C2n^/'.M'nH\A=k). 
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This leads to : 



k=0 ^ ' 



> AcAn:|A|=A: 



cxn'^ ■^\{\A\ =fc)+ I ^ C2n 



< 



fc=0 



A; 



E 7' 

V AcA„:|A|=fc 



< Cr?'^ max 



' ' J fc=0 ^ ^ \AcA„:|A|=A: 



(1^1 = A;) + ^nm = ^) 

C(l-4| = k) 



< Cn^'^ max (7-^) 

ae[A,l] 

Cn^*^ max TT 

,11 -"--^ 



ae[A,l 



n \lx 



Hence, 



n 



'logEi/f (t'^^-^)') < ( n^'^ max V log(a7^ + (1 - a)) ) + o(l) 



as n — )• 00. 



Let us call a™^^ G [A, 1] C (0, 1] the (random) parameter which maximizes the expression above. Then, 
taking the expectation with respect to the environment, and the limit n — )• 00 we get 



lim n^^EelogEff [Y^-^A < hm n^'^E, 



j;iog(ar^(7.-l) + l) 



(24) 



For Bernoulli pinning law, we have to study the function : 



n 



log(a(e/e* - 1) + 1)) + log(a(e/e* - 1) + 1), 



where A^ „ = #{x G A„ : e^; = e} and Nn = i^{x £ An : Sx = e} . It is easy to check that 



(q) < and 5„,(0) = 0. 



^] e+(=f)e-e*) -^0 almost surely 



Moreover, the maximum is attained for a = (g_g^2p(;^_p) \\n: 
as n — )• 00 by the law of large numbers. Hence there exists almost surely no sufficiently large such that 
for n > no, a™^^ = A. Noticing that for a S M"*", gn{ct) < log(A(e/e* — 1) + 1) < 00 almost surely, 
uniformly in n > no, we deduce that 

gn{a) piog(l + A(e/e^-l))) + (l-p)log(l + A(e/e^-l))) 

= Eelog(A(ex/e* - 1) + 1) a.s. and in L\ 

Observe that the r.h.s. is independent of n, so that we can pass to the limit n — )• 00 and obtain 



^ < Eelog(A(e,/e*-l) + l) 
< logEe(A(e,/e*- 1) + 1) 
= log(l) = 0. 
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Let us now do a Taylor expansion of the upperbound (24) around (7 = 0, recalling that cr^ = Yaieiex) 
and that A = C-{d)e*. We obtain: 

p{e)-rie) < plog(l + A(e/e*-l))) + (l-p)log(l + A(e/e*-l))) 
= -c_ {d)K{p)a + 0{a f as o" ^ 



where K{p) = y + P y-\J —f - — y j > for p G [0, 1]. This concludes the proof for d > 3. 

□ 

Case d = 2. The case d = 2 requires some slight modification as i/^ can be stochastically estimated 
by Bernoulli product measures only in a weaker sense. Consequently, we cannot use the same argument 



as in (23). 



By [5l Theorem 2.4, (2.13)] there exists c > such that for all e* sufficiently small, and any set B C A„ 

we have uf{AnB = 0) < (1 - A)l^l where A = A(e*) := ceV log(e'^)| € (0, 1]. 

Then, 



iy^\\A\=k) < u^* {\A\ < k) = uf^* {3B C An ■■\B\=n'^-k and An B = 0) 

k 



< Yl <*(-4ni? = 0)< ("^ )(i-Ar 



BcA„:\B\=n''-k 

We will show that there exists A = A (A) > such that. 



"^^^ (1 - A)"'-'^ < B^d^l = k), for k < [An^ 

Indeed, let k = an'^, we must find > and A(q) such that : 

1_A< AT^(1-A)=:/„(A) (25) 
for all a < oq. It is easy to check that /a(A) has a unique maximum in [0, 1] at A = a. Now, observing 

a 

that fa{oi) = ai-" (1 — a) — )• 1 as a — )• 0, we deduce that for all A > there exists ao(A) G (0, 1] such 
that fa{ci) = i^^'^Aelo 1] faW 1 — A, so that (25) is satisfied. Note that this is sufficient to prove 



/^(e) < /'^(e) for any environment e such that e* > : we have (|.4.| = A;) < (|^| = A;) with 
A > 0, for A; < [Xn'^\ and we can proceed analogously as for d > 3. 

But we need to control how oq varies as a function of A to obtain an estimate on the gap f^{e) — /'^(e). 
As the value A = ce*/y^\ log e*\ is known only for small values of e*, we will only study the case A < 1, 
and find an approximation of the solution to the equation /^(a) = a"/^~"(l — a) = 1 — A around 
a = 0. Developping the left hand side, we get: 

—a + a log a + 0{a^) = —A 

One can check that the behavior of the solution of — a + a log a = — A close to a = is the following : 



is concave on [0,1] (i.e. ^fia*{X)) < 0) and such that ;^/(a*(A))|A=o = —C- Moreover, for C < 1 



A natural guess for the solution of (25) is then a* = — Indeed, /(a*(A)) is a function of A which 
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the solution to the equation /(a* (A)) = 1 — A tends to 1 as C tends to 0. We deduce the following 
result : 

For all Ao G (0, 1) there exists some 6 > sufficiently small such that fa{ot) with a(A) = — satisfies 

(|25]) for all A G [0, Ao]. °^ 

We thus obtain the following upper bound, with A = A(A, 5) = ~]^^ a-nd some fixed e* > sufficiently 

small such that A = A(e*) = ce* / y^l log(e*)| < 1 : 

rie) - r (e) < plog(l + A(e^ 6, c)ie/e* - 1)) + (1 - p) log(l + A(e^ 5, c)(e/e* - 1)) 

+ 0(cj2) as 0- ^ 0, 



c5K{p)a ^ 2^ 



Y^l log e*| log [ce*/^/\\oge 



with K{p) = -^J ^^jp +P (v'S^ ~ \l t\\aX the term log \ ce* / \J\ log e*| j is negative. This 
finishes the proof for d = 2. □ 



3 Attractive/Repulsive environment 

3.1 Introduction 

3.1.1 Description of the results 

In this section, we study a model in which the values of the potential can be negative, so the field 
is penalized if it enters a (small) strip of width 2a around 0. As before, we denote by := 
the set of configurations of the model. Let A d l/. The finite volume Gibbs measure in A for the 
discrete Gaussian Free Field with random "strip'-potential, and boundary condition, is the probability 
measure on O defined by : 



/^A,a('^<^) = exp f -/3HA{ip) + /3 ^ e^.l[<^^g[_a_„]] j n ^'^^ n ^o{dipy). (26) 

A, a \ xgA / xeA yeA" 



where a, /3 £ M"*", and H\{ip) is given by (13). The non-homogenous environment is written e := {ex)x<^A- 
We consider again the law of e given by a Bernoulli product measure. Namely, independently on each 
site a; G A : 

e with probability p , , 

e with probability (1 — p), ^ 

for some fixed triplet {e,e,p). This time e and e can be negative, i.e. e < e G M and p G (0, 1). We 
write again e* and o"^ the average and variance of Cx, respectively (c.f. (|15|)). 

The temperature parameter enters again only in a trivial way. If we replace the field {(px)xeA by 
{^/~^4)x)xeA, o by V^d, and (ex)xeA by {Pex)xeA we have transformed the model to temperature 
parameter (3 = 1. In the sequel, we will therefore work with /3 = 1. 

We define the quenched (resp. annealed) free energy per site in An := {0, ...,n — 1}*^ by (e) (resp. 
/A„(e)): 

flie) := n-^lo, , /^(e) := n-'^log (^^) • 

\ An / \ An / 



We ffist prove in Theorem 3.1 that (/^(e))^ converges to a non-random quantity, denoted by f^{e), 
as A t Z"^ provided |(9A|/|A1 — t- 0. Again, the same proof shows that (/^(e))A converges in M (towards 
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/*^(e)) as A t Z'' provided |5A|/|A| 0. 

The annealed model is the homogenous one with pinning i{e) given by 

£{e) = log(pexp(e) + {I - p) exp(e))) 



(which is different from e* in the (5-pinning case). In Lemma 3.4 we prove that regardless of e both 
the quenched and annealed free energies are non-negative. This motivates the following notions. We 
introduce the quenched (resp. annealed) critical lines, which are delimiting the region where /''(e) = 
(resp. /^(e) = 0) from the region f^{e) > (resp. /'^(e) > 0). We are interested in describing the 
behaviour of these quantities in the phase diagram described by the plane (c, e*) = (y^Var(e), E(e)). 
We hence introduce the so-called quenched and annealed critical lines: 

e*{a) := sup{e* G M : /''(e) = P{a, e^) = 0} and e^^(cj) := sup{e^ G M : r(e) = P{a, e^) = 0} 

Knowing the behavior of the homogenous model for positive pinning |5], we easily deduce that the 
annealed critical line is given by the equation 



;i -p)e V i-f -^pe VP 



£{e) := log(p exp(e) + {1 - p) exp(e))) = ^ e*^(cT) = - log 
Note that /''(e) < f^{e) implies that e*{a) > e*^(fT). We present three results concerning the study of 



the quenched free energy. First of all, in Theorem 3.8 we prove that /''(e) < f^{e) whenever £{e) > 0. 



Secondly in Theorem 3.6 we prove that if the environment is inside {{e,e,p) : /''(e) > 0}, then the 



density of pinned sites in a box of size n is lower bounded by a positive number as n — )■ oo, while in 



the interior of {{e,e,p) : /''(e) = 0}, this density tends to 0. Finally in Theorem 3.9 we show that the 
quenched critical line lies strictly below the axis e* = in the neighborhood of o" = for all d > 2. Our 
result shows in particular that there exists a non trivial region where e* < 0, o" > and /''(e) > 0, 
i.e. where the field is localized though it is repulsed on average by the environment. Unfortunately, 
we don't have any estimate on the behaviour of e*(cr) — e*^{a). Moreover, our theorems do not say 
anything about the existence or absence of smoothing of the transition due to the presence of disorder. 

3.1.2 Extension of pinned sites representation for inhomogenous "strip'-potential 



Following the same lines as Section 2.1.2 for the attractive model, we can perform a high temperature 
expansion of the pinning term : 



^ exp I -l3HA{f) + /3 ^ e^l[^,6[-a,a]] j H ^'^^^ 11 ^o(^'^y'> (2^) 

A, a V x6A / xeA j/eA'= 

l^e^p{-(3HM) J] ((e^^^ - l)l[^,6[-a,a]] + l) n ^'^^ n 



^A,a(^'/') 



y 

A, a xeA zeA yeA': 

^ (t^ ( fie. t\ ^a^V^x e [-a, a] ,Vx G A)\ i w ^ /l^ 

= Mi ^ A^a (^'^ I V'^ ^ '^^ ^ ^) 

AcA VxeA ^A'a / 

V ' 

Observe that the weights can be negative if < 0. For this reason we will only use the measure t'^ ^ 
corresponding to a homogenous pinning potential of strength e > 0. 

For the homogenous case with e > 0, it is known [5J that v^^^ is strong FKG in the sense of [lUj . 
Moreover, it can be compared to Bernoulli product measures in the same fashion as in the (5-pinning 
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case, namely (|5| is still true in this setting. These will be very useful futher. However, let us stress 
that unlike (5-pinning case this measure is not a marginal of fx^^^ giving the "pinned sites". We will 
define it in a next step. 

For some reasons that will be explained later, we will need another representation 



^ exp I -(3Ha{(p) + /3 X] e^llv'^eha.a]] j Yl ^'^^ 11 ^o(^'Py) (^9) 

A, a \ xeA / xeA yeA'^ 



7^ 

A, a xsA xeA S/eA- 

2^ Mi e'^ ^i;o /^A (dv? I ^ = : l^^xl < a}) 

AcA \ \xeA J ^A,a / 



In particular, we obtain a formula for the following ratio of the partition functions : 



,76.0 



7^(e),0 

^A,a AcAxeA 

where ^(e) stands for a homogenous environment for which Cx = for all x £ A, and 7 = {'yx)x£A, 
jx '■= el^'^^'^~^^^'>\ Observe that, contrary to (28), we do not obtain a mixture of Gaussian measures in 



(29). 



3.2 Existence of the quenched free energy 

Theorem 3.1. For d > 2, let f^ie) := | A|-i log(Z^'"/Z"'") be the quenched free energy per site 
A (S Z"^. Then, the limit 



Pie) := lim /^(e) 

AtZ'* 



exists almost surely and in L^, and does not depend on the sequence A t Z'^ provided |M|/|A| 0. 
Moreover, f^{e) is non-random, i.e. 

rie)=E,ine)) a.s. 



Proof. The idea of the proof is the same as for Theorem 2.2 so we only write the part that must be 
adapted, namely the upper-bound on the height of the field under /Xn'a. We keep the same notations 
as previously. 

Lemma 3.2. There exist Ci, C2, C3 > such that for all x £ An and for all n sufficiently large, 

l^faiW-\ >T)< Cie-^^(^-^3logn)Vlogn_ 



Proof. We first use (29) to get 



l^tA\^^\>T) = J2'^nAA).f,'/j\^x\>T\A = {x:\ipx\<a}). 

AcA 

As it was mentioned before ^jf'{dip \ A = {x : \^Px\ ^ a}) is not Gaussian making the proof considerably 
harder. We will use the results listed in survey |12| . We are going to show that 

\ A = {x : \ipx\ < a}) I G A„, (px>a), 
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where -< denotes the stochastic domination. To this end we define, for k > : 

f^fcj-^ ^) iH^x + a)%^^^[_a,o]]+k{^x-a)H[^^^[Q^a]] for X G A, 

^ [H^x - a)'^l[^^>„] + k{ip^ + a)'^l[^^<_„] for x £ 

Further we define also measures fJ-fi^p) oc ex.p{—Hn{^) — ^^x&A^ ^ii^x)) for i G {1,2}. One checks 
that condition |12| (B.3)] is fulfilled hence -< /^l- Passing with k — >• +00 we get the aforementioned 
domination and as a consequence 

I^A^ifx > T\ A = {y. \ipy\ < a}) < fJ-°/^J<fx > T | Vy G A„, ipy > a). 

For d > 3 we use |12| Theorem 3.1], which in our notation implies 

An := fJ^A^i^x I Vy G A„, (py > a) < C^/logn, 

for some C > 0. By the Brascamp-Lieb inequality [12| Section B.2] we obtain 

fJ'A^iy^x > T \\/y e An, ify > a) < - A„ > T - ^„ | Vy G A„, > a) 

< eM-Ci{T-Anf/logn) 

for another constant Ci > 0. 

The case d = 2 is slightly more involved. Our first aim is to prove that An < Clogn, for some constant 
C > 0. Let us assume on contrary that An > Clogn. By |12| (B.14)] we have fij^^{ipf) < logn + A'^. 
Using the Paley-Zygmund inequality we get 

> AJ2) > \ . > 1/8, 

for n large enough. Now choosing C large enough we get a contradiction of f3l Theorem 4]. The rest 
of the reasoning follows now the same line as before. 

□ 

Again, exactly the same proof suffices to show the following 

Theorem 3.3. For d>2, let /^(e) := |A|-Mog(EeZ^'°/Z°'°) he the annealed free energy per site in 
A (E Z"'. Then, the limit 

r(e) := lim /X(e) 

exists in R and does not depend on the sequence A f Z'^ provided |5A|/|A| — )■ 0. 

□ 

3.3 Bounds on the quenched free energy 
3.3.1 Posit ivity of the quenched free energy 
Fact 3.4. For d>2, 

r(e) := lim n-^^log M^f >0 
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Remark 3.5. This proves also that /'^(e) > since by the Jensen inequality f^{e) < /'^(e). 
Proof. For < e < e (and any d> 1) the proof is easy, since by the Jensen inequahty 

= ""'i^g/^A'! ( ( E «-i[b.i<«] ) ) ^ i™o""'^A! ( E «-i[iv'.i<«] ) ^ 

V / / \a;eA„ / 

The non-trivial part of the claim concerns e < (no restriction on e which is by convention bigger 
than e). We will prove that /^(e) > erj for all rj > 0, with different strategies for d = 2 and d>3. 

Proof for d = 2. In this case, 

0>0 /ur I I / 1 ^ d\ ^ /^a! (H^ ■\^x\<a}) 1 ^ 0,0. 1 I ^ X ^ 2a 



a;eA„ ' V n 

The last inequality comes from the fact that ipx is a Gaussian variable under /x^'^ with variance 

'^n = )"a°(</'i) ^ Clogn ^ oo as n ^ oo. The right hand side is then less than 1/2 for n sufficiently 
large, hence 

\ \xeA„ / / \ V xeAn / / 

from which we deduce /'^(e) > eij. Finally we set rj ^ 0. 

Proof for d>3. In dimension 3 and above, the variance of the free field is bounded, so the previous 

argument does not work. But for a cost of surface order wc can shift the boundary condition to a 
constant M > 0, and from this height the probability of having a density of sites where \(px\ < a which 
is bigger than rj is less than a half for M sufficiently large. Indeed, recalling that we write A*^'^ for the 
measure with boundary condition M, for some C > we have 



\ \ x&A„ / / \ \ x&A„ / 



Now, 

/^A„_i [^{^ ■ \'Px\ < a} > r?(n - 1) J < ^^jin^ " ^(^ _ i)^ E ^'An-^\^^\ ^ °) 

= , ^ ,w y thf {iPx>M-a) 

^ \ 0,0 / 2\ 

- {M - a)^r]{n - l)d 2^ ^An-i^*^^^ 

< ^F}'^^ ^o < 1/2 for M sufficiently large. 

r]{M — ay 



Hence, 



which implies as before that /'^(e) > ery for all r/ > 0. □ 
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3.3.2 Density of pinned points 



Theorem 3.6. Let e be the environment as described in (21), the following holds : 



1. If, additionally, e is such that f'^{e) > 0. Then there exists some c > such that for almost 
every realisation of e we have 



lim Jn-'^ Vl 
\ ieA„ 



[</3iG[-a,a]] > C 



2. Let Z := {e €z {e, e}^'' : /''(e) =0}. IfeEZ, then for any c > and for almost every realisation 
of e we have 

■ " (30) 



\ ieA„ / 



where Z denotes the interior of the set Z. 



Proof. The proof of the first claim goes similar lines of the proof of Fact 2.6 and hence is skipped. 
For the second claim note that (e) is a function of e and e. Let us write 



Pn-=n ^ l[e,=e]l[^,e[-a,a]] 

xeA„ 



It is easy to prove that 



3.1 



that it 



^/a„(^) e,0 /- n , ^^/a„(^) e,0 /- 

—Q^=l^A„JPn), and g-2 = Var^„,a(Pn) > 0- 

Hence (/^ (e))n is a sequence of convex functions in the variable e. We know by Theorem 

o 

converges to /''(e). Therefore, for e & Z we have 

r e,o , dflie) g/q(e) 

The same reasoning holds for p := X^xeA^ -''-[ei:=e]l[(/Pa:e[-a,a]]- The second claim is now a straight- 
forward consequence of the Chebychev inequality. □ 

Remark 3.7. The question of the behaviour on the boundary of the 0-phase is intimately connected 
with the question of the degree of the phase transition which stays unresolved for the moment. 

3.3.3 Strict inequality between quenched and annealed free energies 

Theorem 3.8. Take {p,e,e) such that p £ (0,1) and e <e (no restriction on the sign of e). Let e 



be an i.i.d. Bernoulli environment of parameters {p,e,e) as described in (27). We recall the notations 
e* = pe+ {1 — p)e and cr^ = (P'p + e^(l — p) — e*^. 

(a) Let d> 3. Then, 

/"(e) < /-(e). 
Moreover, there exists C3{d,a,e*) > such that, 

/-(e) - /''(e) > C3(d,a,e*)cj2 + 0(a3) as a ^ 0. 
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(b) Let d = 2. There exists e > uniform in p,e,e such that if e^^'^') — 1 < e, then 

r(e) < r(e). 

Moreover, under the same hypothesis on i{e), there exist C2(a,e*) > such that, 

f''{e)-f'^{e)>C2{a,e*)a^ + 0{a^) as a ^ 0. 



Proof. The same idea as the proof for the attractive 5-pinning case can be As in Section 2.3.2 , we have 



7e,0 



r(e) - r(e) = hm n-'^Eelog = hm log(P^(^)(7-^)) 



Z 



A„,a 



with £{e) = log(pexp(e) + (1 — p) exp(e))) is the annealed environment corresponding to e, and = 

The following simple estimation is crucial for our proof. We are able to upperbound P^'-'^^d^l = k) 
using i/^'-'^^d^l = k) and consequently, using results of [5j, by a Bernoulli measure 

AcA„ 

AcA,i,\A\<k AcA,i,\A\<k 



with A := C-(a)(e^^''^ — 1). Let us notice that the above much resembles (23). The proof follows further 



unchanged up to (24) where we have to analyse the function 

+ (1 - a)) + log (a • e'^-^(^) + (1 - a)) . 



5n(a) := 



log [a-e 



We can easily check that < and gn(0) = 0. Moreover, the maximum is attained for 

~ (ee^e^Xi-pjp ^~ ^ ~ ~^ ~^ *-* almost surely as n — ?• oo by the law of large numbers. The 

same reasoning holds then, and we end up with 

gn{a) Ee log(A(e"--^('') - 1) + 1) a.s. and in L\ 

leading to the upper bound on the difference of free energies 

r(e)-r(e) < Eelog(A(e^^-^(-)-l) + l) 
< logEe(A(e^--^(^)-l) + l) 
= log(l) = 



(31) 



using the strict convexity of the logarithm and the definition of ^(e). The Taylor expansion of the 
upperbound (31) around o" = gives the following estimate. (We recall that o"^ = Vare(ea;), and that 
A = c_(a)(e^(e)-1).) 



r(e)-r(e) > [l[e''* -l^\_{af]a' + 0{a') as a ^ 
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This concludes the proof for d > 3. 

The generahzation for d = 2 follows the same lines, with A = C-{a){e^^''^ - 1)/ V| log(e^('=) - 1)|. We 
obtain : 

This concludes the proof for d = 2. □ 
3.4 Attraction by a repulsive on average environment 

Theorem 3.9. Take (p,e, e) such that p G (0,1) and e < < e. Let e he an i.i.d. Bernoulli envi- 



ronment of parameters {p,e,e) as described in (27). We recall the notations e* = pe + (1 — p)e and 
0-2 = e2p + 6^(1 -p) - e*^. 

For d > 2, the quenched critical line is located in the quadrant {{a, e*) : o" > 0, e* < 0}. 

More precisely, for d> 3, there exists some C,C' > depending on d, a,p only and e G (0, 1) such that 
for any environment e which fulfills e > 0, — e < e < and 

fe*>r^r% ford = 2 

) log{-e) 



e* > -C- for d>3, 



we have f^{e) > 0. 



Remark 3.10. 1. Both of this bounds can be seen on Figures ^ in the plane (a,e*). Moreover, 
the bound for d > 3 can be rewritten as e* > —C"{d, a,p) ■ a^. 



2. Jensen's inequality gives us an upper bound on C,C'. Indeed, as f^{e) > f^{e), if f^{e) = 
then /''(e) = 0. In particular, we must have —C < ^^e.c\e=o ^ '^^'^ result gives thus an 
upper-bound on the behavior of the quenched critical line near cj = 0. 



Figure 2: Phase diagram of the model. The red curve is the annealed critical line; the two other curves 
are the bounds we get on the position of the quenched critical line for d > 3 (blue curve) and d = 2 
(green curve). They are both valid for small ex. (Note that the green curve could be above the blue 
one, depending on the constants C and C", which we do not control. ) 
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3.4.1 Proof for (i> 3 



Proof. The idea is to tilt the measure such that the field ip is shifted up of an amount h on the sites x 
for which ex < 0- In this way the shift of the field follows the environment. For some technical reasons, 
we need to work with the measure with boundary condition a, so perform two changes of measure (first 
one changing boundary condition and the second one to following the environment). Let h > (to be 
fixed later). 



/1(e) log /^A^,e,/^ ( . O^a " /o.^ ^P [^.6[-a,a]] 



where {'^x)xi^Kn ^ii^er is distributed as {ipx + ^l[e:r<o])^g^ under /U^'^. More formally, intro- 

ducing T^^h ■■ {{Vx)xeAj ^ {'P^ + ^l[e.<0])^.6A„' ^^^^^ ^^A^e,h /^A,! ° ^eT^ 
Using Jensen's inequality, we get 



/ 0,a J 0,0 

/^A„,e,h exp 2^ ea[^,e[-a,a]] + log — + log 



d/i 



An,e,h 



A 0,a 



\ 0,a 



1 0,a 1 0,0 

E ^^H^.^\-aA\ + log + l°g 7% 



xeA„ 



v 



(1) 



(2) / 



As Z 



0,0 



^A eh~ (which follows by change of variables in the Gaussian integral), the first term can 
be written as ^ 



{x,y}nA„^0 



where ipx '■= '^x + ^l[ej,<o]- Hence, using the definition of //^'^ ^ ^, 



^-Vte,.((l)) 



n 0,a 



{x,?/}nA„^0 

x^y 



Ad 



Y (l[e.<0] - 1k<0])' 



{a;,)/}nA„^0 
xr^y 



\ 



J 



The second term contains only boundary contribution of order n'^ ^. Indeed, 

^ A 



(2) = l2a o^ldAnl + log ^ > 2a ^ - a^iaA^I - Cri 



xedAn 

where we used the same methods as in Section 



An 



xedAn 



d-1 



2.2 



to estimate the ratio Z^"'/Z^^. Hence, 



x£dAn XGdAn 
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We get 



/A„(e) > n ^ e^fij^y^^ £ [-a, a]) 



x£A„ 



4^d ^ -[^-<o] 

{x,?;}nA„^0 

0,a 



Now we use the fact that the marginal laws of all (/^j;, x G A„ under [i^ are Gaussian variables centered 



at a, i.e. ^px ~ A/'(a, cr,^) where o"^ = [ipx) < VarJ^'*((/7a;) < c((i) < oo for d > 3. Therefore, 



An 



/^a!('/'^ ^ [-«' "]) - ^a!(¥'^ + ^ ^ «]) = /^a1'(¥'^ ^ [-2a, 0]) - ii'^JyHix G [-2a - /i, -/i]) 



0,a 



0,0, 



0,0, 







C , 

'— /i J — 2a—h 

h as n — )• oo, 



-2a 



(32) 



for c((i) » h. In particular we will use that : 

fJ'A^^x e [-a, a]) - f^^y^Px + he [-a, a]) > Ci{d, a) • h, 
for some Ci((i, a) > 0. 

/„^(e) > n-" {t^T^. G [-a, a]) - a)/il[e,<o] , 



x-eA„ 



4d 



E (l[e.<0] -l[e,<0])'-Cn-l 



{x,y}nA„7^0 



Observe that fij^"'{ipx G [—a, a]) = iijf{ipx G [— 2a, 0]) > fi^{ipx G [— 2a, 0]) > C2{d,a) for some 
C2{d,a) > 0. By taking the expectation with respect to the environment, using the bounded conver- 
gence theorem and Theorem 3.1 we get : 



Pie) = lim Ee/^^ (e) > e*C2id, a) - hCi{d, a)pe - h^p{l - p) /4 



(33) 



We may optimize over h as the left hand side does not depend on it. Doing this one checks that 
/'1(e) > as soon as 



Ciid,a) 



p 



--: —K{d, a,p) ■ 



C2{d,a) 1-p 

This gives the implicit equation in terms of the variance o"^ of Cx ■ 
1 



1 



4a 



'^~2K^''\I~^2 ^ Kp^ 



K{l-p)a^ 
P 



+ 0{a' 



The annealed critical curve as well as this bound are drawn on Figure [2j We recall that (32) is valid 
under assumption that h is small. The maximum of (33) is realized at /imax = ' thus it is 

enough to assume that e is small. 

□ 
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3.4.2 Proof ford = 2 



In the case d = 2, the variance of the Gaussian free field diverges with the size of the box, so we cannot 
use the previous estimates. To circumvent this problem we introduce the so-called massive free field. 
Let m > 0, 



z 



0,i 

An,m 



exp 



\ 



8d 



2 2 

m 



\ 



{x,y}r\Anj^0 



xeAn 



XGAn x£dAn 

(34) 

Known facts about this model can be found in [71 Section 3.3]. In particular, the random walk 
representation introducted for the massless GFF is still true, but for a random walk Yj that is killed 

2 

with rate ^{m) = j^^, namely at each time i, if the walker has not already been killed, the process 
is killed with probability ^(m), where the killing is independent of the walk. 

Lemma 3.11. Let d = 2. Then, 

1. There exists some Ci > such that for n large enough, m > small enough and all x E A„, 

^JiZ,^{vl)<Cl\\og{m)\. 

2. There exists some C2 > such that for n large enough and m > small enough, we have 



,0,0 



n loe 



Z 



0,0 

An,0 



> — C2?7T-^| log(m) 



Proof. These bounds are rather standard. We give here the main steps of the proofs with some 
references. For the first claim, we use the random walk representation [7] to write 

00 00 

l^Z.mivl) = E ^-(^^ = ^ ' A H > ^) = ^(1 - i{m)fP^{Yi = X , TA„ > ^) 
£=0 1=0 

where ta^ is the first exit time of A„ and H is the killing time of the random walk Yt. Hence, 



(35) 



e=i 



where Xi is a simple random walk (without killing). The projections of Xi onto the two coordinate 
axis are two independant 1— dimensional random walks Xj and Xf, then by Stirling formula. 



PoiX^e = 0) = {Po{Xl, = 0)f 



2^ \ 2~2e 



1 



(1+0(1)) as £-^00 



The asymptotics of (35) for small m gives the desired upperbound. 

To prove the second claim we use the representation of the partition function described in [4, p. 542] (it 
applies to the massive GFF with an obvious modification). We denote by P the coupling of a random 
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walk Xn and a killed random walk Yn such that = Xn up to its killing time 



|An| Mog 



Z 



0,0 
A„,0 



Z 



0,0 



(1 oo 1 

2 E E ^ (^-(^2£ = X, TA„ > 21) - P,{Y^, = X, rA„ A H > 2£) 

^ 2 E 27 = 0, TA„ > 2^) - Po(i"2£ = 0, rA„ A H > 2^) 

oo 

= 2E^^o(^2^ = 0'^A„ >2£,b^ <2^) 

^ oo 



2 ^ 2£ 



(36) 



Using the same estimate as in (35), the asymptotics of (36) for small m gives the desired upperbound. 

□ 

The idea is to tilt the measure, as in the proof for d > 3, first to work with the massive measure, and 
second to follow the environment such that the field (p is shifted up of an amount h on the sites x for 
which Cx < 0. For some technical reason, we need to work with the measure with boundary condition 
o, so we perform three changes of measure (first one for changing boundary condition, a second one 
for adding mass, and a third one for following the environment). 
Let h > and m > to be fixed later. 



0,0 



/aJs) = n log^^' 



exp ^ exl[ 



n log MA„,m,e,h 2^ exl[^.e[-a,a\\ + Tl)^ + TlM^ + log 



where i^x)xeA„ under fJ^Al,m,e,h distributed as (93^. + /il[e^<o])^g^^ under ; more formally, 

introducing T^^h ■ {{fx)xeAj ^ + ^l[e.<0])^6A„' ^^^"^ /^Alm,e,h as /^Alm^^eTh- Using 
Jensen's inequality we get 



(e) > n-V^': 



J 0,0 J 0,s 1, 0,s 

2^ exl[0,s[-a,a]]+log-^+log--^+log 

d^A„ d^A„,m d^A„,m,e,h 



zGA„ 



V 



(1) 



(2) 



(3) 



As in the proof for d > 3, we have 



By Lemma 



3.11 



we have 



(2) = log 



Z' 



yO.S 

0,S - 

A„,m 



«-V^':™,e,.((l)) > -Cn~\ 



^0,3 ^0,0 ^0,0 



^n:^ > -Cn - Can^m^l logm|, and then 

^A„,m ^A„,0 '^An.O 



0,s 
A„ 



+ m y. V^x ^ — C'n — m \ logm| + m y, 



X&An 



xeA,i 



28 



hence, 

™~VA':,™,e,J(2)) > -Cn'^ - C2rr?\ logm|. 

Finally, noticing that -^A^meh ~ '^^K m 0^^^ perform a change of variables in the Gaussian integral), 
we can compute the third term. 

(3) = -^ ^ {H^x-H^yf -{(px-(pyf -rn^^{}Px-sf -{(p^-sf 

where (px '■= + ^l[e^<o]- Now we will use the fact that the marginal laws of all ip^, x S A„ 
under Ij^^^ ^ are Gaussian variables, i.e. tpx ~ M{iJL^,a^) where fi^ ^ s except for x close to the 
boundary of the box. Indeed, by the random walk representation of the mean, there is C > such 
that l/iA'.mCv'x) - s| < C(l + ^7^2)-||x-aA„||^ Moreover, {a^f = Var°'^ < Ci| logm|. Using the 

definition of /i^'^ ^ ^ ^, and computing the terms as in the proof for ci > 3, 

{x,y}nAnyt0 x£A„ 
x^y 

We get, for n large enough and m small enough 

/^Je) > n-'^ E ex/^Atm(^^ ^ [-a, a]) ^ ^ (l[e.<o] - l[e,<o])^ 

2 2 2 ^ ^ / 2 1 1 

—n~ m h y, l[e^<o] — C m \ logm| — Cn~ , 

xeAn 

for some C and C > 0. Note that for 0{n'^) sites x, we have 

f^Al,miv^x e [-a, a]) - fJ'Al^rni'Px + ^ ^ I""' "]) ^ ^1,^2 («) ' ^ as n ^ oo, (37) 

for h « a < s = Ci|logm|, and a = Ci|logm|. Above ^s,a stands for the p.d.f. of the above 
Gaussian distribution with mean s and variance a^. In particular, for a positive fraction of x (close to 
1) and m sufficiently small, we have the upper bound : 

f^Almifx e [-a, a]) - ^A'lm(¥'x + ^ e [-a, a]) > y^^^ • h, 
for some Ci(a) > 0. Now we can compute : 

f^Je) > ^ ex{fi1;U^x G [-a,a]) - |^^^^l[e.<0]) 

x&A„ ' ^ ' 

^ Yl (l[e.<o] - l[e^<o])^ - n-^m^/i^ E l[e.<o] " C'm'^\ logm| - Cn"^ 

{x,j/}nA„7^0 a;eA„ 

Observe that /^a'I m('/'a; ^ a]) > 2o • $^^2(— o) = |^g^^| uniformly in x G An- Let us take the 
expectation with respect to the environment, and use the bounded convergence theorem and Theorem 
3.1[ we get : 



r(e) = hm Ee/^ (e) > - h ■ ^f^^ - h^m^p - ■ - C'm^\ log m\. {38) 

" ' ~- " ' logm| I logm| 
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Our aim now is to show that the right hand side can be positive even when e* is negative. In the 
above expression h, m are free parameters which we may vary. However, we have to remember that 
both h and m need to be small enough, which makes standard optimisation analysis cumbersome. We 
are going to show that there exists C > and Bq < such that for any e E (cq, 0) and 



:= C, 



log(-e) 



there exist small h and m such that the rhs of (38) is positive. Notice that the result will imply that 

2 

for any e* > Cy^^gp^ the free energy is positive. Let us choose the value of h which maximizes (38 1 
for fixed m, i.e. 



h 



Ci{a)e 



and for m let us take 



(2m2p+ (1 -p)/2)|log(m)|' 
^2 - -k/{\ogkf, where k := -e*Ci{a)/C'. 



(39) 



m~ = -fc/(logfcJ", wnere := -e"Ui[a)/U' . (40) 

One can verify that with the above choice of parameters both h and m are as small as we want. Let 
us first put (39) into the rhs of (pSj) and obtain 



C'rn^l logm| 



C'm^l logm|. 



logm| (4m2p + (1 — p))(logm)2 
For k and consequently m small enough we have 

|iogm| 2(1 — pj(logmj^ 
Further let us multiply both sides by (logm)^ and insert (40). 

/'i(e)(logm)2 > -e*Ci{a)logm + ^^j-"^^^- +C'm\logmf 

2(1 -p) 

= -"^^(log A; - 31og(| logfcD) + ^^^^ - C'^^^,{logk- 31og(| logA;|))= 



.C'i(a) , , , s CiiaVpe"^ ^^(a) / *n , * 

-e*^^logi-e*Ciia)/C') + ^^^^^ + e'^^ + oie*) as |e* 







From the last claim it is straightforward to conclude existence of C (sufficiently small) and Bq with the 
properties described above. □ 



Acknowledgements. We would like to warmly acknowledge Yvan Velenik for introducing us to the 
topic, and for many usefull discussions. L.C. thanks Fabio Toninelli for his advice. 
L.C. was partially supported by the Swiss National Foundation. P.M. was supported by a Sciex 
Fellowship grant no. 10.044. 



30 



References 



[1] K. S. Alexander. The effect of disorder on polymer depinning transitions. Comm. Math. Phys., 
279(1):117-146, 2008. 

[2] K. S. Alexander and V. Sidoravicius. Pinning of polymers and interfaces by random potentials. 
Ann. Appl. Probab., 16(2):636-669, 2006. 

[3] E. Bolthausen, J.-D. Deuschel, and G. Giacomin. Entropic repulsion and the maximum of the 
two-dimensional harmonic crystal. Ann. Probab., 29(4): 1670-1692, 2001. 

[4] E. Bolthausen and D. loffe. Harmonic crystal on the wall: a microscopic approach. Comm. Math. 
Phys., 187(3) :523-566, 1997. 

[5] E. Bolthausen and Y. Velenik. Critical behavior of the massless free field at the depinning tran- 
sition. Comm. Math. Phys., 223(l):161-203, 2001. 

[6] P. Caputo and Y. Velenik. A note on wetting transition for gradient fields. Stochastic Process. 
Appl., 87(1): 107-1 13, 2000. 

[7] A. Dembo and T. Punaki. Lectures on probability theory and statistics, volume 1869 of Lecture 
Notes in Mathematics. Springer- Verlag, Berlin, 2005. Lectures from the 33rd Probability Summer 
School held in Saint-Flour, July 6-23, 2003, Edited by Jean Picard. 

[8] J.-D. Deuschel and G. Giacomin. Entropic repulsion for the free field: pathwise characterization 
in d> 3. Comm. Math. Phys., 206(2) :447-462, 1999. 

[9] J.-D. Deuschel and Y. Velenik. Non-Gaussian surface pinned by a weak potential. Probab. Theory 
Related Fields, 116(3):359-377, 2000. 

[10] C. M. Fortuin, P. W. Kasteleyn, and J. Ginibre. Correlation inequalities on some partially ordered 
sets. Comm. Math. Phys., 22:89-103, 1971. 

[11] S. Galluccio and R. Graber. Depinning transition of a directed polymer by a periodic potential: 
A d-dimensional solution. Phys. Rev. E, 53:R5584-R5587, Jun 1996. 

[12] G. Giacomin. Aspects of statistical mechanics of random surfaces. Notes of the lectures given at 
IHP, 2001. 

[13] G. Giacomin. Random polymer models. Imperial College Press, London, 2007. 

[14] G. Giacomin, H. Lacoin, and F. Toninelli. Marginal relevance of disorder for pinning models. 
Comm. Pure Appl. Math., 63(2):233-265, 2010. 

[15] G. Giacomin, H. Lacoin, and F. L. Toninelli. Disorder relevance at marginality and critical point 
shift. Ann. Inst. Henri Poincare Probab. Stat, 47(1):148-175, 2011. 

[16] G. Giacomin and F. L. Toninelli. Smoothing effect of quenched disorder on polymer depinning 
transitions. Comm. Math. Phys., 266(1):1-16, 2006. 

[17] D. loffe and Y. Velenik. A note on the decay of correlations under (5-pinning. Probab. Theory 
Related Fields, 116(3):379-389, 2000. 



31 



[18] E. Janvresse, T. de la Rue, and Y. Velenik. Pinning by a sparse potential. Stochastic Process. 
AppL, 115(8):1323-1331, 2005. 



[19] C. Kiilske and E. Orlandi. A simple fluctuation lower bound for a disordered massless random 
continuous spin model in D = 2. Electron. Comm. Probab., 11:200-205 (electronic), 2006. 

[20] C. Kiilske and E. Orlandi. Continuous interfaces with disorder: even strong pinning is too weak 
in two dimensions. Stochastic Process. AppL, 118(11):1973-1981, 2008. 

[21] F. Spitzer. Principles of random walks. Springer- Verlag, New York, second edition, 1976. Graduate 

Texts in Mathematics, Vol. 34. 

[22] A. C. D. van Enter and C. Kiilske. Nonexistence of random gradient Gibbs measures in continuous 
interface models in d = 2. Ann. Appl. Probab., 18(1):109-119, 2008. 

[23] Y. Velenik. Localization and delocalization of random interfaces. Probab. Surv., 3:112-169 (elec- 
tronic), 2006. 

[24] Y. Velenik. Le modele d'Ising. Lecture notes, 2009. 



32 



